The fate of vacancy-induced supersolidity in 4 He 
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The supersolid state of matter, exhibiting non-dissipative flow in solids, has been elusive for thirty 
five years. The recent discovery of a non-classical moment of inertia in solid 4 He by Kim and Chan 
provided the first experimental evidence, although the interpretation in terms of supersolidity of the 
ideal crystal phase remains subject to debate. Using quantum Monte Carlo methods we investigate 
the long-standing question of vacancy-induced superflow and find that vacancies in a 4 He crystal 
phase separate instead of forming a supersolid. On the other hand, non-equilibrium vacancies 
relaxing on defects of poly-crystalline samples could provide an explanation for the experimental 
observations. 
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The observation of a non-classical moment of inertia in 
solid 4 He by Kim and Chan (KC) has provided the first 
experimental evidence of a possible supersolid phase of 
matter [J, , which is characterized by crystalline order 
and frictionless flow. Early theories of supersolidity were 
based on the assumption that the low-temperature He 
crystal may be incommensurate (the number of atoms 
is not an integer multiple of that of lattice sites). As 
a consequence of their quantum behavior, point defects 
such as vacancies and interstitials can Bose condense at 
low temperature, giving rise to superflow. 

In particular, the Andreev-Lifshitz-Chester (ALC) sce- 
nario 0, 0] assumes that the gain in kinetic energy by 
delocalizing the vacancy can overcome the potential en- 
ergy cost of creating it in a perfect crystal, such that a 
dilute gas of highly mobile vacancies can be stabilized. 
The energetics of this scenario is illustrated in the lower 
curve in Fig. ^ Similar supersolid ground states with a 
low density of strongly correlated vacancies were recently 
put forward by Dai, Ma and Zhang (DMZ) [j| and also 
by Anderson, Brinkman and Huse (ABH) [fj. Their sce- 
narios are possible even when single-vacancy excitations 
in the perfect crystal are gapped, as indicated in Fig. ^ 
It was also recently suggested that the onset of superso- 
lidity leads to anomalies in the elastic moduli 0. 

Although the initial experimental observation by KC 
has been confirmed by other groups H, 0, [ljj, the in- 
terpretation in terms of superflow in the crystal phase 
is becoming increasingly questionable. Repeated cycles 
of annealing make the supersolid signal weaker to van- 
ishing Measurements of pressure-driven flow have 
yielded a null result in hexagonal-clo se-p acked (hep) 4 He 
|lll IT^ I . Experimental evidence 0, ll-l| points toward a 
commensurate hep ground state in 4 He, as the measured 
vacancy activation energy is approximately 10 — 15 K, 
and appears to rule out thermally activated vacancies at 
the characteristic temperatures below 0.2 K of the KC 



experiment. X-ray measurements [lH ] on hep 4 He put 
the tightest upper bound on the vacancy concentration 
n v (T -> 0), disfavoring the ALC and DMZ/ ABH sce- 
narios, although ABH debate the interpretation of these 
measurements @ . On the theoretical side, two of us have 
proven that a superfluid crystal must generically be in- 
commensurate Il5j . a view which is supported by recent 
numerical simulations showing that ideal hep crystals of 
4 He are insulating pL6l 1171 fia | . 

The physics of vacancies in solid Helium is an impor- 
tant open problem, even aside from its relevance to the 
experiment of KC, and has been studied for a long time. 
The majority of microscopic calculations have been vari- 
ational |l9l I20I l2l| , and focused on the properties of a 
single vacancy. Various estimates of the vacancy activa- 
tion energy |19ll20| . including the one obtained by Ceper- 
ley and Bernu using Path Integral Monte Carlo simula- 
tions fTfj| . are in quantitative agreement with experiment. 
However, Galli and Reatto questioned the reliability of 
all previous numerical calculations, including their own 
ones [2l]], by raising the importance of finite size effects. 

Here we present results of a Quantum Monte Carlo 
study of vacancies in hep 4 He at low temperature (0.2 K 
< T < 2 K). For single crystals in the thermodynamic 
limit, we find large vacancy and interstitial activation 
energies Ay and At at all densities, ruling out the ALC 
scenario. We then proceed with the discussion on finite 
vacancy concentrations to see if activated vacancies form 
a (meta) stable gas (upper dashed line in Fig.QJ, or if they 
support the DMZ/ABH-type state with a small vacancy 
concentration n\. It turns out that the only minima 
in the grand-canonical (free) energy are at zero vacancy 
concentration, corresponding to a commensurate solid, 
and at high vacancy concentration such that the system 
is in a liquid phase. In Fig. ^ this scenario is called "un- 
stable" . 

We employ the worm algorithm |22j , a grand-canonical 
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FIG. 1: (Color online) A sketch of the different scenarios for 
vacancy states in the ground state of solid Helium at the melt- 
ing density. Supersolidity can occur at finite vacancy densi- 
ties nv, according to the Andreev-Lifshitz-Chester (ALC) or 
Dai-Ma-Zhang(DMZ) /Anderson-Brinkman-Huse ( ABH) sce- 
narios. ALC assume an energy gain already for a single va- 
cancy (negative initial slope), while DMZ/ABH postulate a 
ground state with a low density of strongly correlated vacan- 
cies, without making any prediction for the single vacancy 
behavior. For the scenarios with gapped vacancies (positive 
slope Ay), a (meta)stable gas of out-of-equilibrium vacan- 
cies will be formed if the energy curve bends up (the upper 
dashed line) due to repulsion between the vacancies. The 
curve labeled 'unstable' represents the situation supported in 
this paper: An attractive vacancy gas is unstable and phase 
separates into the commensurate insulating crystal (with den- 
sity ns) and the liquid (with density «l < Us). 



Path Integral Monte Carlo technique formulated in 
the configuration space of the single-particle Matsubara 
Green's function. In our simulations, we use a three- 
dimensional cell, designed to fit a perfect finite crystal of 
N 4 He atoms, with periodic boundary conditions in all 
directions. Our microscopic model of 4 He is standard, 
based on the accepted Aziz pair potential J^. Single- 
particle properties can be deduced from the Matsubara 
Green function 

G(t,t',t) = (T^r^ViO)}}, (1) 

which can easily be sampled in the worm algorithm. Here 
(...) stands for the thermal expectation value, T is the 
time-ordering operator, — j3/2 < t < (3/2, ?/^(r',0) and 
"0(r, t) are Bose particle creation and annihilation oper- 
ators, respectively. For r > (r < 0) one is computing 
the Green function for an interstitial atom (vacancy). 

In Fig.|21we show the spatially averaged Green function 
G(k = 0, r) for the hep 4 He crystal at melting density 
n = 0.0287 A~ 3 and for a low temperature T=0.2 K. Re- 
sults are shown for a system of ./V=800 atoms. Since va- 
cancies and interstitials are gapped by a finite activation 
energy, the low temperature Green function G(k = 0, r) 



FIG. 2: (Color online) Single-particle Green function G(k = 
0,t) computed for hep 4 He at the melting density n o =0.0287 
A -3 and T — 0.2 K. Symbols refer to numerical data, solid 
lines are fits to the long-time exponential decay. The given 
numerical values are the interstitial (Ai = 22.8 ± 0.7 K) and 
the vacancy (Ay = 13.0 ± 0.5 K) activation energies, inferred 
from the slopes of G. The inset shows the vacancy-interstitial 
gap -Egap = Ai + Av for different system sizes. 

decays asymptotically as 

G(k = 0,r) ~e-l T l A , 

with activation energies A = Ay (Ai) for t < (r > 0), 
which are measured relative to the chemical potential fi 
of the system. In order to determine the chemical poten- 
tial at the melting- freezing curve, fj, , we have computed 
by simulation the equilibrium density of the liquid, n([x), 
at low T. We located the experimentally known freezing 
point n(// ) = 0.02599 A~ 3 with high accuracy. Our es- 
timate is fi = 0.06 ± 0.04 K. Upon fitting the observed 
exponential decays of G(0,t), we obtain Av=13.0(5) K, 
and Ai=22.8(7) K. The value of A v found here is only 
slightly lower than results from other numerical calcula- 
tions, whereas that of Ai is a factor of two smaller than 
the one reported in the only comparable study [16j. 

We did not find any appreciable size dependence of ei- 
ther Ay or Ai, as seen in the inset of Fig. [2J Since the 
results are also temperature independent for T < 1.5 K, 
our data yield reliable ground state estimates in the ther- 
modynamic limit. We can thus confirm that single vacan- 
cies and interstitials in a perfect 4 He crystal have huge 
activation barriers at the low temperatures of the KC 
experiment. This finding is consistent with the absence 
of thermally activated vacancies at low T, as reported in 
X-ray experiments [l3| . 

The same Green function analysis was extended to 
higher densities. From Fig.[3|we infer that both Ay and 
Ai increase with density, or equivalently, pressure. This 
can be understood intuitively, as the cost for creating 
a defect in a denser system is higher. Large activation 
energies thus prevent the ALC scenario in solid Helium 
at any density. It seems reasonable to extend this con- 
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FIG. 3: (Color online) Chemical potential thresholds for cre- 
ating a single interstitial (upper curve, /ii) and vacancy (lower 
curve, /xy) at different densities. The middle curve, /ihcp is 
the chemical potential of the hep solid deduced from the com- 
pressibility data and the value of Mhcp = M° a t melt- 
ing density. The gaps A V (i) for activating a vacancy (or in- 
terstitial) are given by the difference in chemical potentials, 
Mhcp - A*v (or fii - ^hc P )- 



elusion to all presently known, naturally occurring solids 
since they have stronger interatomic potentials. 

We proceed to the physics of multi-vacancy states to 
check whether the highly correlated DMZ/ABH ground 
state in Fig. ^ can be formed, despite the large activa- 
tion energy Ay- We also explore the possibility of a 
(meta)stable gas of out-of-equilibrium vacancies. The 
thermodynamic stability of a uniform dilute gas of va- 
cancies can be addressed by performing simulations in 
the grand-canonical (fi — V — T) ensemble. We have per- 
formed simulations at low temperature T = 0.2 K with 
N = 800 and 2016 atoms in the crystal under different 
pressures up to 65 bars. We manually created vacancies 
by explicitly removing a number of atoms. Various initial 
configurations with randomly located, remote and clus- 
tered vacancy positions were considered, with vacancy 
concentrations ranging from half a percent up to six per- 
cent. 

Our consistent observation is the formation of a com- 
mensurate hep crystal, as the grand-canonical simulation 
fills the vacancy sites with atoms. This behavior is seen 
not only when the chemical potential fi is set equal to the 
experimental chemical potential /i = Mhcp, but also when 
it is shifted several K below /ihep, which can be explained 
by the large activation energy Ay- When the chemical 
potential /j, is lowered below that of the liquid but is still 
kept above the single- vacancy threshold ( /iv < M < Mo), 
a runaway in vacancy number for all finite initial concen- 
trations occurs, leading to a liquid phase. 

We conclude that a uniform dilute gas of vacancies 
is thermodynamically unstable against separation into 
a vacancy-rich and a perfectly crystalline vacancy-free 
phase, which does not melt: The scenario corresponding 
to the curve labeled "unstable" in Fig. 2] is realized in 
4 He, no DMZ/ABH- type state is found. The instability 
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FIG. 4: (Color online) The vacancy-vacancy correlation func- 
tion u(r) as a function of the distance r between the vacan- 
cies shows that three vacancies easily cluster and form a tight 
bound state. The inset shows a typical snapshot of a layer of 
atomic positions averaged over the time interval [0, f3] (filled 
black circles), where we see that the three vacancies (trian- 
gles) have a tendency to cluster in layers. 



occurs not only at a finite vacancy concentration, but 
also for a few vacancies in a large crystal. 

In order to confirm this conclusion, we performed sim- 
ulations in the canonical (N — V — T) ensemble, keep- 
ing the number of particles fixed. The instability trans- 
lates into phase separation, with the vacancies forming 
clusters inside a perfect crystal. Indeed, our simulations 
produce the spontaneous formation of vacancy clusters 
for temperatures up to 1 K, regardless of pressure and 
the initial configuration. Phase separation is also sup- 
ported by energy calculations, since the energy difference 
(E(ny) — E(0))/ny is not monotonically increasing with 
vacancy concentration ny. 

Inspection of the simulations revealed that already 
three vacancies cluster easily and form a tight bound 
state, as shown in the inset of Fig. 0] The vacancy- 
vacancy correlation function is the quantity most sen- 
sitive to the sign of their effective interaction. We have 
computed it by first averaging atomic world lines over 
an imaginary-time scale of the order of 0.2 K _1 , which is 
much longer than the zero-point atomic motion but much 
shorter than the vacancy-atom exchange time. Next, we 
compare the obtained atomic positions with the lattice 
points in an hep crystal. Using a similar procedure as 
in Ref. [l^, we repeatedly remove the closest particle- 
lattice point pairs from the list, until we end up with a 
small number of unmatched lattice points which we de- 
fine as the vacancy positions. The decay of the vacancy- 
vacancy correlation function v{r) in Fig.0]shows that the 
vacancy-vacancy interaction is attractive and that three 
vacancies form a tight bound state. A similar study for 
two vacancies also yielded attractive correlations but was 
less conclusive regarding a bound state; if it exists, it is 
rather shallow with a binding energy of about IK. 

The observed attractive interaction between vacancies 
- leading to phase separation - arises from the elastic 
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interactions. These are long ranged, decaying as 1/r 3 , 
and are attractive in certain directions (in the basal plane 
and along the c-axis) due to dominant quadrupolar terms 
[2(j. A similar effect was previously su gges ted for 
bound states of substitutional 3 He atoms [23,|28|, and it 
is also reminiscent of a proposal for oxygen vacancies in 
high-temperature superconductors [2£j • The strong effec- 
tive short-range attraction is due to the local minimum 
in the pair potential. 

Given the instability of the vacancy gas, the only 
remaining possibilities accounting for the experimental 
data of KC are inhomogeneous scenarios, due to a less 
than ideal sample quality 0, |3(j ■ Metastable vacancies 
can play a pivotal role in this case. When a large va- 
cancy cluster is formed, the system can macroscopically 
lower its energy by deforming the cluster into a disloca- 
tion loop. We found, by simulating artificially created 
edge dislocations in the hep crystal, that these can be 
either insulating or sustain superflow depending on the 
annealing procedure for the initial setup. In particular, in 
the presence of a large stress field close to the dislocation 
core, the core first melts locally into a liquid, which then 
forms a highly disordered superfluid (an undercooled su- 
perglass [l7|h 

Other scenarios include the migration of vacancies to- 
wards the grain boundaries and boundary edges, along 
which superflow might be possible: imagine, as a 
schematic example, an insulating grain boundary, where 
the potential produced by the crystal bulk forces the 
atoms at the interface to form an insulating Mott state. 
When vacancies relax on the grain boundary, they effec- 



tively dope the Mott insulator and may make it super- 
fluid — a phenomenon well documented in lattice models. 

In conclusion, we have shown that vacancies are unsta- 
ble in a 4 He crystal. They form clusters and the system 
phase separates into a vacancy-rich phase and aperfect, 
insulating crystal. Recent experimental results [8j are in 
perfect agreement with the fact that the ground state 
of solid 4 He is thus never a supersolid, but a commen- 
surate crystal. The experimental observation of a non- 
classical moment of inertia in solid Helium is due to non- 
equilibrium quantum behavior, such as superflow along 
defects or the formation of a superglass around crystal- 
lographic defects. High resolution structural analysis of 
defects in the 4 He samples and large scale simulations of 
elementary defects will be important in completing our 
understanding. 
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